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ABSTRACT 

In the standard scenario of the fireball model of gamma-ray bursts(GRBs), 
the huge initial energy release produces a relativistic blast wave expanding 
into the external medium and a reverse shock moving into and heating the 
fireball eject a. We calculate the high energy gamma-ray emission due to inverse 
Compton scattering of the synchrotron photons from relativistic electrons in 
the reverse shock. Under the favorable values of the physical parameters of the 
GRBs and the interstellar medium, our result shows that during the prompt 
phase, this emission dominates over the component from the forward shock at 
high energy gamma-ray bands. This mechanism can excellently account for the 
observations of the prompt high energy gamma-rays detected by EGRET, such 
as from GRB930131. 



Subject headings: gamma rays: bursts — radiation mechanisms: non-thermal 



1. Introduction 



Gamma-ray bursts (GRBs) are widely believed to be caused by the dissipation of the 
kinetic energy of a relativistically expanding fireball with a Lorentz factor of i] ~ 10 2 -10 3 
(see e.g. Piran 1999 for a review). After producing the prompt GRB, the shell impacts on 
the surrounding gas, driving an ultra-relativistic shock into the ambient medium, which 
approaches a self-similar behavior (Blandford & McKee 1976) after a short transition phase. 
This forward shock continuously heats fresh gas and accelerates electrons, producing a 
long-term afterglow. The initial interaction between the ejecta with the surrounding gas 
also leads to a reverse shock which propagates into and decelerates the ejecta. This shock, 
heating up the shell's matter and accelerating its electrons, operates only once. Once the 
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reverse shock crosses the ejecta, the ejecta expands and cools adiabatically. Transition of 
the fireball to the self-similar behavior occurs on a time scale comparable to the reverse 
shock crossing time of the ejecta. Thus, emission from the fireball ejecta is suppressed after 
the transition to the self-similar expansion. 

A prompt 9-th magnitude, optical flash has been recently detected (Akerlof et al. 
1999) accompanying GRB990123. The most natural explanation of this flash is synchrotron 
emission from a reverse shock propagating into the fireball ejecta after it interacts with 
the surrounding gas (Sari & Piran 1999a; Meszaros & Rees 1999). The reverse shock, at 
its peak, contains energy which is comparable to that of the GRB itself and the forward 
shock, but has much lower temperature than that of the forward shock. So it radiates 
at considerably lower frequencies. Since, during the prompt phase, the number of heated 
electrons in the reverse shock is 10 2 — 10 3 ) times larger than in the forward shock and 
the seed photon source is stronger, we expect a strong synchrotron self-Compton emission 
flux at high energy gamma-ray bands in the reverse shock. Here we present an analytic 
calculation of this synchrotron self-Compton emission and consider its implications for 
the prompt high energy gamma-ray emission detected from some GRBs by the Energetic 
Gamma Ray Experiment Telescope (EGRET) and Compton Telescope (COMPTEL) on 
board Compton Gamma- Ray Observatory (CGRO). Early discussions on the synchrotron 
self-Compton process in the reverse shock include the works of, e.g., Meszaros et al. (1993) 
and Sari et al. (1996), both of which, however, are quite different from our subject in many 
fundamental aspects. 

2. The synchrotron self-Compton emission from the reverse shock 



Transition of the fireball ejecta to the self-similar expansion occurs on time scale At tran 
comparable to the longer of the two time scales set by the initial conditions: the (observed) 
GRB duration AtcRB and the (observed) deceleration time scale td ec of the fireball ejecta, 
at which the self-similar Lorentz factor (r) equals the initial ejecta Lorentz factor 77, i.e. 
A-Uran = max [ A£g.r£ , tdec] (Sari & Piran 1999b; Waxman & Draine 2000). 

The deceleration time tdec is in fact the observer time at which the heated interstellar 
medium (ISM) energy is comparable to the initial energy E, i.e. 



2.1. 



Reverse shock model 



"('dec 




(1) 
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where 

r dec = ( , f 2 ) 1/3 = 3.7 x lO^cm^fno 1 / 3 ^ 3 (2) 

is the deceleration radius, E = 10 53 i?53erg, n = lnocm -3 , 77 = 300i] 300 , and z is the redshift 
of the GRB source. Since even long duration GRBs have durations of order of lOsec, which 
is comparable to tdec, we therefore take At imn ~ tdec for the following calculations. 

During the transition, the unshocked fireball shell propagates at the original 
expansion Lorentz factor rj and the Lorentz factor of the plasma shocked by the 
reverse shock in the rest frame of the unshocked ejecta is 7 rs ~ i]/T. According to 
Sari & Piran (1995), the reverse shock Lorentz factor 7 rs depends on the parameter 

£ = (//A) 1 /2 r? -4/3 = 2E yS n - l/3 A -l/2 % -4/ 3j where l = ( _J_)l/3 the Sedoy length &nd 

A = cAt GRB is the shell width. If £ 3> 1, the reverse shock is Newtonian and 7 rs ~ 1, while 
if £ 1, the reverse shock is relativistic and 7 rs 1. Since the shell is likely to spread, 
which adjusts £ > 1 to £ ~ 1 (Sari & Piran 1995), the reverse shock we considered here will 
be always mildly relativistic (7 rs ~ 1) for typical parameter values of the fireball shock. 

The distributions of the electrons in both the forward and reverse shocks are assumed 
to be a power-law of index p (^(7) oc 7~ p ), with the Lorentz factor of the random motion 
of a typical electron in the shell rest frame being 7 m = ^^elsh, where £ e is the fraction 
of thermal energy carried by electrons (Waxman 1997a,b ). Here 7^ = V for the forward 
shock while 7^ = 7 rs for the reverse shock. Assuming that is the fraction of the thermal 
energy carried by the magnetic field, we have the peak frequency of the reverse shock 
emission 

•C = = 1.6 x io«H Z (iy^)^j M „; /2 (-^-), 0) 

1 + z 2irm e c 0.3 (J.Ol 1 + z 

where B' = 12G(^-) 1 / 2 77 3 oo^o^ 2 , is the magnetic field in the comving frame. Here £ e and t; B 
are the values relevant for the reverse shock and have been assumed to be not different from 
the forward shock (Wang, Dai & Lu 2000). The numerical values we have used are those 
characteristic of the forward shock (Waxman 1997a; Freedman & Waxman 2000; Wijers & 
Galama 1999; Granot et al. 1999). As this peak frequency we estimated is slightly higher 
than the optical band and it decreases with time after the deceleration time td ec , we expect 
that the optical light curve should posses a short initial rise phase, and then decays with 
time, which is consistent with the light curve of the prompt optical flash from GRB990123. 
Waxman & Draine (2000) have also derived the peak flux for the reverse shock 

r- = 1 . 5Jy ^ 5( -^)^)./ 2 , Ji „-v.^ ( _|_w_ r3/4> (4) 
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for a flat universe with zero cosmological constant and H = 65/i 65 Km s~ 1 Mpc~ 1 . This flux 
is consistent with the observations of GRB990123 within a factor of two, whose peak flux 
is about Uy (Akerlof et al. 1999). The peak frequency of synchrotron emission from the 
forward shock is larger than that of the reverse shock by a factor ~ (T 2 /r]) 2 , and is, at the 
deceleration time tdec of the ejecta, 

i/£ = 1.4 x l0 20 Hz(^) 2 (i^) 1 / 2 % 4 00 no 1/2 (^-), (5) 

v o.3 ; v o.or /300 v i + z' y ' 

since at this time the self-similar Lorentz factor r ~ i]. 

The cooling break frequencies , at which the cooling time scale of electrons equals that 
of the expansion dynamic time, are given by 

_ 10 17 Hz £ B 3/2 4 3/2 t dec 2 l + z 
Uc ~ (Y rs + ir [ 0M } r]30 ° n ° 40sec j ( 2 j W 

for the reverse shock, and 

j 10 17 Hz £ B ^_ 3 /2 4 -3/2 ( tdec \-2/ j^ + Z \ f 7 \ 

vi = (r^Ti? W } 77300710 W ( ^~ } (7) 

for the forward shock, where Y rs and Yf s are the Compton parameters for the reverse shock 
and forward shock, respectively (Panaitescu & Kumar 2000). The Compton parameter 
Y, expressing the cooling rate of electrons due to inverse Compton effect, is defined 
as Y = |r e / 7 2 iV e (7)(i7, whereiV e (7) is the normalized electron distribution and r e is 
the optical thickness to electron scattering. According to Panaitescu & Kumar (2000), 
Yf s = Ky^II^ + 1 — 1) ~ 2, since the forward shock is in the fast cooling regime, while 
Yrs ~ (7m) P ~ 1 (7c S ) 3 ~ Pr e S ~ 3 for the reverse shock, where r r e s is the optical thickness of 
the shocked ejecta to electron scattering, which will be calculated in the next subsection 
and 7 C is the Lorentz factor where the total energy emitted at the time t^ ec is comparable 
to the electron's energy: lc = = fPg Vo^)" WK The 

synchrotron spectrum of the forward shock is thus described as // s = f^iy I v{ s Y 1 ^ 2 for 
> v > v{ s and // s = f^(^ jv{ s Y x l 2 (y '/v^)~ p l 2 for v > while for the reverse 
shock, ft' = f™(v/v™)- P -^ for i/" > v > and ft = fZ KV^T^KT 1 fOT 
v > v r c s . 



2.2. The synchrotron self-Compton emission from the reverse shock 



The synchrotron self-Compton emission from the reverse shock can be calculated by 
using the above equations for the synchrotron spectrum and the Thomson optical depth of 
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the shocked fireball material, which is 

r?=^- = 8-4 x 10-«£ 5 f nrvS^r" 2 , (8) 
^r dec t dec 

where N is the total number of electrons in the fireball and <jt is the Thomson scattering 
cross section. This inverse Compton luminosity for the reverse shock is larger than the 
synchrotron luminosity by a factor Lie / L sy ~ Y rs ~ 3 , dominating the synchrotron cooling 
rate of the electrons, and radiates at much higher energy band than the synchrotron 
emission. We need consider only the first-order inverse Compton and neglect the higher 
order inverse Compton process, because a once-scattered synchrotron photon has energy of 
order (/iz/ m j n )com7f > Tn e c 2 in the rest frame of the second scattering electron with Lorentz 
factor being 7 e and thus the Thomson limit no longer applies. Since the reverse shock is in 
the slow cooling regime, the up-scattered spectrum peaks at 

^ W = 2(7^) 2 ^ = 1.0 x 10 21 Hz(||) 4 (^ r )V^/ 2 Of)\^-l (9) 

and the cooling break frequency, at the deceleration time t dec , is 

v rs,C = 2{l r c s fv rs = ^ x l 22 Hz(i^)^/ 2 % 2 / 3 n - 13 / 6 E 5 - 3 4/3 ( T ^), (10) 

where we have substituted the expression of td ec (Eq.l)into it. At this deceleration time of 
the ejecta, the peak flux of the inverse Compton component is 

fZL° = < s f- = 2.6 x 10- 8 erg cm^s^MeV-^^no 7 / 6 A 3 (^) 1/2 /4( / • 

(11) 

For the simplicity of the analysis, we use the broken-power law approximation to the inverse 
Compton component, though the presence of logarithmic terms makes some corrections 
to the spectrum for v > v r ^ lc (Sari & Esin 2000), i.e. f™> IC = f r r ^ IC \v I v'^ 10 )- 2 ^ for 

v rs,IC > y > v rs,IC and frsJC = frs,IC ^rsJC j^ICy^- ^ J y rs,ICy\ for v > u rs,IC _ p or 

the sake of comparison, we give the flux of the inverse Compton component at three 
representative frequencies: 

f rs ' IC (e = 5MeV) = 2.6 x 10 _8 erg cm-V^eV^^ 3 ; 
frs,W( £ = 10 0MeV) = 2.7 x 10" 9 erg cnrV^eV- 1 ^ 3 ; 
f rs ' IC \e = lGeV) = 1.5 x 10~ 10 erg cm- 2 s- 1 MeV- 1 E 5 2 3 /3 , (12) 



where we have only retained the relations with the shock energy E 53 and replaced other 
shock parameters and the redshift with the typical values. 
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After the reverse shock has passed through the ejecta, the ejecta cools adiabatically. 
Sari & Piran (1999a) assumed that the ejecta follows the Blandford-McKee (1976) 
self-similar solution, in which a given fluid element evolves with a bulk Lorentz factor of 
7 ~ i? -7 / 2 , and successfully explained the decaying light curve of the optical flash and 
the radio flare behavior of GRB990123 (Kulkarni et al. 1999). Below we will give the 
decaying light curve of the synchrotron self-Compton component of the reverse shock after 
it has passed through the ejecta (t obs > td ec ). Since oc t~ b J^ 48 ~ t~ bs and r r e s oc t ^J 2 , 
the peak flux of the inverse Compton spectral component f™ ,IC = fm T l s ^ Kb!" 2 ■ If 
the observed frequency locates between the two break frequencies {y r c s ^ IC < v < v r ^ IC ), 
then f™> IC = f£JC(_v iu )-(p-iy2 m According to Sari & Piran (1999a), 7m oc t~ b f /4S and 

vZ oc t~ ob T\ so v^ lc = HrVZ oc t~ b f /w ~ t- b 2 s , thus f^ IC oc t'^ oc t£ for p = 2.5. 
On the other hand, if the observed frequency is above v r c s > IC ', the flux drops exponentially 
with time since all the electrons above the corresponding energy cool and no fresh electrons 
are accelerated once the reverse shock has crossed the ejecta shell. Therefore, in general, we 
expect to see a rapidly decaying high energy flux from the reverse shock, which constitutes 
a unique characteristic distinguished from other models (see the next subsection) suggested 
for the observed high energy gamma-ray emission from some GRBs. Measurements of the 
time dependence of the high energy gamma-ray flux with the planned Gamma-ray Large 
Area Space Telescope (GLAST) mission will test this synchrotron self-Compton scenario. 



2.3. Detectablity of the prompt high energy gamma-rays and the case of 

GRB930131 

EGRET has detected 'prompt emission above 30MeV from seven bright GRBs triggered 
by BATSE (Catelli et al. 1998), among which 1.2GeV and 3.4GeV photons have been 
detected from GRB930131 (Sommer et al. 1994; Ryan et al. 1994) and GRB940217 
(Hurley et al. 1994), respectively. GRB940217 also displays delayed high energy gamma-ray 
emission, about 90 minutes after the initial trigger. There are a handful models suggested 
to explain the delayed and prompt GeV emission. For examples, Katz (1994) suggested 
that the impact of the fireballs with dense clouds (n > 10 9 cm~ 3 ) could yield high energy 
gamma-ray emission via it decay; Waxman & Coppi (1996) proposed that the cascading of 
ultra high energy (> 10 19 eV) cosmic rays (accelerated in cosmological GRBs) off infra-red 
and cosmic microwave background fields can produce delayed GeV-TeV emission due to 
the deflections by the intergalactic magnetic field; Dermer et al. (2000) argued that the 
synchrotron self-compton emission from the forward shock expanding into a rather dense 
circumburst medium (n ~ 100cm -3 ) may be responsible for the prompt and delayed high 
energy gamma-ray emission; Vietri (1997) and Totani (1998) suggested that the synchrotron 



-7- 



emission of the protons may be responsible for the GeV-TeV emission. Here we suggest 
that the synchrotron self-Compton emission from the reverse shock could explain both the 
flux level and the spectrum of the high energy gamma-rays detected by EGERT. As an 
example, we discuss the case of GRB930131. The high energy gamma-ray emission from 
GRB930131 extends for about 25s. Totally sixteen gamma rays above 30MeV were detected, 
including two with energy about lGeV, and the power-law fit of their photon spectrum 
is given by g ~ 7.4 x l(r 6 photons (cm s MeV)- 1 (5/147MeV)- 2 - 07±a36 (Sommer et al. 
1994). In the last subsection, we have found that the inverse Compton emission intensity 
from the reverse shock is f^ IC ~ 2.7 x l(T 9 erg cm^s^MeV-^/lOOMeV)- - 75 ^ 3 
for e < hv r c s ' IC and f r v ^ IC ~ 2.7 x l(T 9 erg cnrt^MeV-^/lOO MeV)" 1 - 25 ^ 3 
for e > hv r c s > IC with p = 2.5, that is, the differential photon flux is 
| ~ 6.0 x l(T 6 photons (cm s MeV)- 1 (5/147MeV)- 2 - 25 E 5 2 3 /3 . Thus, if the fireball 
shock energy E ~ 1.5 x 10 53 erg and other parameters such as £ e , rj, z and the number 
density n of the surrounding medium take the above representative values, then both the 
flux level and the spectrum agree well with the observations. However, since this emission 
process lasts only a time comparable to the crossing time of the reverse shock, it cannot 
explain the delayed high energy gamma-ray emission, as observed from GRB940217, which 
may, for example, be due to deflections of ultra high energy cosmic rays accelerated in the 
GRB fireballs by the intergalactic field (see Waxman & Coppi 1996). 



2.4. Comparision with the high energy gamma-ray emission flux from the 

forward shock 

For completeness, we also compute the high energy gamma-ray emission flux from the 
synchrotron and self-Compton process in the forward shock, and compare them with that 
from the reverse shock. The peak flux of the forward shock is 

/£ = 8mJy/. 2 5 ( v ^-^ i ) 2 (ig T )V 2 E 53n y 2 (13) 

(Waxman 1997b; Wijers & Galama 1999), and the two break frequencies v£ and 
are given by Eq.(5) and Eq.(7). Please note that here the flux peaks at u^ s . For typical 
parameter values of GRBs, during the prompt phase the forward shock is in the fast cooling 
regime, and we get the high energy flux at those three representative frequencies: 

f fs (e = 5MeV) = 1.2 x l(T 8 erg cm" 2 s" 1 MeV" 1 J B 53 ; 
f fs (e = lOOMeV) = 2.5 x l(T 10 erg cnrt^MeV -1 ^; 

ff s (e = lGeV) = 1.4 x 10" n erg cm- 2 s- 1 MeV" 1 J B 53 . (14) 
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We can see that the flux at ~ 5MeV is comparable to the inverse Compton component of 
the reverse shock, while at higher energies (> lOOMeV) it is about one order of magnitude 
lower. As in the reverse shock, we also need consider the synchrotron self-Compton emission 
process in the forward shock, since the typical energy of the synchrotron photon in the 
rest frame of the scattering electrons is of the order of {hvl s ) com r ) e ^.m e c 2 , which is still in 
the Thomson limit range. At the deceleration time (t b s ~ td ec ), the characteristic break 
frequencies of the inverse Component of the forward shock are 

= 2( 7 £« = 8.0 x ^°M^n^) 1/2 v!oonon^ z ) (15) 

and 

u fsjc = 2{li s fv f s = 5 x 10 22 Hz(^)-^^^ (16) 
The peak flux of this inverse Compton component is 

ft JC = r! s fl s ~ VTnrfl: = 3 x l(T 13 erg cm^MeV^^E^^-^nf . (17) 

i dec U.U1 

Therefore the inverse Compton flux from the forward shock at e ~ 250MeV band is 
~ 3 x 10~ 13 erg cm~ 2 s _1 MeV _1 , much lower than that from the reverse shock. They are 
comparable only at e ~ lOTeV band. So, we conclude that for the typical parameter values 
of the shock and the surrounding medium, the synchrotron self-Compton emission from 
the reverse shock dominates over the synchrotron and synchrotron self-Compton emissions 
from the forward shock at high energy gamma-ray bands. Our result is different from that 
of Dermer et al. (2000), who argue that the synchrotron self-Compton emission from the 
forward shock may be responsible for the prompt and delayed high energy gamma-ray 
emission. The key point of the difference is that they considered a rather dense circumburst 
medium with number density n ~ 100cm -3 , while we consider a typical interstellar medium 
with n ~ lcm~ 3 . 



3. Summary and Discussion 

The detection of the prompt optical flash and the radio flare from GRB990123 has given 
positive evidence for the reverse shock model. In this paper, we calculated analytically the 
synchrotron self-Compton emission from the reverse shock. Assumed favorable parameter 
values for the fireball shock and the density of the interstellar medium to be £ e = 0.3, 
£b = 0.01, E 53 = 1, p — 2.5, rj = 300 and n = lcm~ 3 , the result shows that during the 
prompt phase, this emission dominates over the emission due to other processes of the 
shocked electrons in both the forward and reverse external shocks at high energy gamma-ray 
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bands. We suggest that this process provides a new explanation for the prompt high 
energy gamma rays detected by EGRET from some bright bursts, such as GRB930131, 
GRB910503 and GRB940217 etc. As an example, we compared our calculated result with 
the observations of GRB930131 and found that both the flux level and spectrum agree quite 
well, provided the initial fireball shock energy is a few times 10 53 erg, which is consistent 
with the fact that all the bursts detected by EGRET are bright bursts in the BATSE range. 
We also derived the decaying light curves of the inverse Compton component and found 
that it decays quite rapidly, regardless whether the observed band locates above or below 
the cooling break frequency of the inverse Compton component. This unique feature can 
act as a test of our model if in the future GLAST is capable of describing the time evolution 
of the high energy gamma-ray flux. Because once the reverse shock has crossed that ejecta, 
all electrons cool by adiabatic expansion and no fresh electrons are accelerated, our model 
can not explain the delayed GeV emission, as observed from GRB940217. 

GRBs are generally believed to be produced by the internal collisions between the 
relativistic shells at a radius r ~ 10 12 — 10 13 cm, which is much smaller than the distance at 
which the external shock takes place. Till now, the mechanism of prompt GRB emission is 
not as well-understood as that of afterglows and we do not know whether the usual GRBs 
and high energy gamma-rays have a common origin . If we here take a naive assumption 
that the internal shock emission extends from BATSE band to GeV bands with f v oc z/~ p / 2 , 
the simple extrapolation gives a flux of f u (e = lOOMeV) ~ 1.0 x lCT 8 erg cm~ 2 s _1 MeV _1 for 
typical GRB fluence of F 7 (20 KeV - 2 MeV) ~ l(r 5 erg cm" 2 . This flux at e = 100'MeV 
is comparable to that from the reverse shock. However, this simple extrapolation may be 
incorrect if there is a sharp cutoff at some high energy band, as shown in the model of 
Ghisellini & Celotti (1999), in which the prompt GRB emission is due to quasi-thermal 
comptonization process in internal shocks and ends with an exponential cutoff at energies 
higher than a few MeV. 

In all our calculation, we considered a fireball shock expanding into a uniform 
interstellar medium with constant density of n ~ lcm~ 3 . The observations of some 
afterglows have shown that some GRBs may occur in a stellar wind environment (Dai & 
Lu 1998; Chevalier & Li 1999), implying the massive star origin of these bursts (Woosley 
1993; Paczyhski 1998). As having been pointed out in Chevalier & Li (2000) and Dai & 
Lu (2000), the peak flux of the reverse shock formed by the initial interaction between the 
fireball ejecta with the stellar wind medium is significantly lower than that in the case of 
the interstellar medium for similar shock parameters . For this reason, we conjecture that 
the high energy gamma-ray emission in this kind of reverse shock is too weak to be detected 
by high energy gamma-ray detectors. 
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